POISSON BRACKETS OF ORTHOGONAL POLYNOMIALS 



MARIA JOSE CANTERQi AND BARRY SIMON^ 



Abstract. For the standard symplectic forms on Jacobi and CMV matrices, 
we compute Poisson brackets of OPRL and OPUC, and relate these to other 
basic Poisson brackets and to Jacobians of basic changes of variable. 



1. Introduction 

It has been known since the discoveries of Flaschka [M] and Moser [32] concerning 
the finite Toda lattice as a completely integrable system that finite Jacobi matrices 
of fixed trace support a natural symplectic form. Explicitly, if 

/6i ai \ 



J = 



V 

on the manifold where aj > and 



(1.1) 



a.N-1 



N 



the nonzero Poisson brackets of the a's and &'s are 
{bk, ak} = -T flfc 



fc = 1, 



.N 



(1.2) 



(1.3) 
(1.4) 



For the basic definitions of symplectic manifolds and Poisson brackets, see Deift [5]. 

The nonzero pairings are neighboring elements in J. As one would guess, there 
are differing conventions, so that sometimes the j in p.3p and (|1.4p are ^ and 
sometimes one takes { , } to be the negative of our choice. 

It is easy to check that X^jLi commutes with all a's and 6's as is necessary 
for this to define a bracket on the manifold where (|1.2p holds. It is also easy to 
see that the form defined by (|1.3p / (|1.4p is nondegenerate and closed so there is an 
underlying symplectic form; see the end of Section [31 
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If one takes 

N N-1 

H = 2TiiJ^) = 2^b^+i^a^j (1.5) 

j=i i=i 

then the Hamiltonian flow, 

bk = {H,bk} = 2ial~al^,) (1.6) 

flfe = {iJ, flfe} = afe(6fe+i - 6fe) (1.7) 

is the Toda flow in Flaschka form [14]. It is completely integrable; indeed, (see 
(frnj) below) 

{Tr(J"),Tr(J™)} = (1.8) 

for all n, m. 

If dp is the spectral measure for J and (1,0,..., 0)*, then 

N 

The E^nd {xj}^^^ are not independent; indeed, by (|1.2p . 

=c ^p,-l (1.10) 

Fundamental to developments are the following Poisson brackets: 

{xj,Xk} = Q {xj.pk} = \[5jkP'j - PjPk] (1-11) 

The analog of {xj^Xk} = for periodic Toda chains is due to Flaschka [14j and 
(jl.lip appeared (implicitly) first in Moser [32 . 

Our original motivation was to understand the analog of this for CMV matrices. 
Along the way, we realized that Poisson brackets of orthogonal polynomials were 
not previously studied, even in the Jacobi case, so we decided to discuss both cases. 

We also felt that the centrality of (|1.11|) . which is partly known to experts, was 
not always so clear, so we also decided to say something about that. Indeed, one 
of our initial goals was to prove the analog of {xi,Xj} = for the periodic case of 
OPUC whose previous proofs [41] Section 11.11] and 34j[35] were involved. Hence 
one application of p. lip we discuss is going from (|l.lip to the periodic case. 

Since ^'^d X^jLi Pj ^-re constant, {xj, p^} must sum to zero over j or k, 

as can be checked in (jl.lip since Y^^=i Pj — ^■ 

We note one curiosity of (|l.lip . namely, {xj,pk} is symmetric in j and k. We 
will eventually also find 

{p^^p^} ^ .liE^ _ y PiMii. + y MhEi^ (1.12) 

^ - ^ I ' ^ nit . nit ' ^ nit , ivt 

•^1 t^k / ■ -^1 -^fn I , t^k -^m 

but will make no use of this complicated formula here. It is needed if one wants to 
find angle variables for the Toda flows (see [2^ ITB] )- 

Simultaneous to our work, Gekhtman-Nenciu |16j were studying Poisson brack- 
ets of Caratheodory functions, which we will see is closely related to our work. 
Indeed, a key to our progress was learning of this work and the earlier paper 
of Faybusovich-Gekhtman [13|. We also mention not unrelated earlier papers of 
Kako-Mugibayashi [26l[27] and Kulish [31]. 
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For the monic orthogonal polynomials on the real line (OPRL), Pn{x), our basic 
result is that ior n — 1,2, . . . , N, 

{Pn{x),Pn{y)} = {Pn-l{x),Pn-l{y)} = (1.13) 
' Pn{x)Pn-l{y) - Pn{y)Pn-l{x) 



2{P„(a;),P„_i(2/)} = 



y 



- Pn-l{x)Pn-l{y) (1.14) 



Notice the occurrence of the Bezoutian (see [5D]) of Pn and Pn-i in the first 
term on the right of (|1.14p . (|1.13p and (|1.14p will follow by a direct induction in n. 

There is a symmetry between P„, P„_i and P„, Qn (where Qn is the second kind 
polynomial (of degree n — 1) for P„), for P/v is a determinant of z — J and Pn~i 
(resp. Qn) are the determinants obtained by removing the last and rightmost row 
and column (resp. first and leftmost row and column). This will immediately lead 
to 



{Pn{x),Pn{y)} = {Qn{x),Qn{y)}=0 
2{P„(x),Q„(y)} - - Pn{y)Qn{x) 



x-y 



(1.15) 

Qnix)Qn{y) (1.16) 



The sign changes because inverting order changes the sign of { , }. All these calcula- 
tions appear in Section [21 In [T3|, several different Poisson brackets on polynomials 
are considered, including one that is essentially equivalent to (|1.16p . 

In Section [21 we derive {xj,Xk} — from {P„(x), P„(y)} = 0, {xj,pk} from 
{Pnix),Qniy)} and {pj,pk} from {Qn{x),Qn{y)} = 0. 

Looked at carefully, if we argued directly for P, Q by induction, we, in essence, 
use coefficient stripping to relate Pn.Qn for {aj}^^^ U {bj}^^^ to Pn-i,Qn-i 



for {aj}^S'2 U {bj}^S'2- The analog for orthogonal polynomials on the unit circle 
(OPUC) is a new coefficient stripping relation of first and second kind paraorthog- 
onal polynomials, which is new. We present this in Section [4l denoting these first 
and second kind OPUC by P/v and Qn- 

In Section [5l we use the natural symplectic form for OPUC: in terms of Verblun- 
sky coefficients, 

{aj,ak} = Q {aj,dk} = -ip'jSjk 
This was introduced by Nenciu-Simon [37] ■ We will find 

{P„(z), P„(?«)} = = {Qn{z), Qn{w)} 



(1.17) 



(1.18) 



iPn{z)Qniw) ~ Pniw)Qn{z)) 



z + w 



w 



Qn{z)Qn{w) + Pn{z)Pn{w) 



(1.19) 



using induction and the coefficient stripping of Section [H P„ , Qn will depend on a 
parameter fi. {ajY^II^ — > {~ctj}j'=o which preserves OPUC Poisson brackets and 
j3 —f3 interchanges P„ and Qn, so the right side of (|1.20|) has to be antisymmetric 
under interchange of P„ and Q„, as it is. 

In Section [6l we derive the analogs of (|1.11|1 . In (|1.11|) . the x's and p's are global 
variables fixed by xi < a;2 < . . . . For OPUC, the eigenvalues are only locally well 
defined, namely, 



N 



(1.20) 
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where 



(1.21) 



The analogs of (jl.lOp are 



N 



N 



c 



(1.22) 



and of prn|) . 



{e,,ek} = o 



(1.23) 



The analog of {9j,6k} = for the periodic case appeared first in [37], and 
{ej,ek} = is from [31 [35]. The formula for {9j,Hk} is due to Killip-Nenciu 

In Sections [THTTl we discuss three applications of the fundamental relations (|l.lip 
and (|1.23[) : to computations of Jacobians of the maps (a,b) {x,p) and (a) 
{0,ii), to the exact solution of the flows generated by G{x) or G{9), and to the 
periodic case. In Sections [Hi and [TUl we discuss the differential equations induced 
on the monic OPRL and OPUC. Sections [T^ and [T31 compute more Poisson brackets 
for OPRL and OPUC, respectively. 

It is a pleasure to thank Rafael Hernandez Heredero and Irina Nenciu for useful 
discussions. M. J. Cantcro would like to thank Tom TombrcUo and Gary Lorden 
for the hospitality of Caltcch. 

2. Poisson Brackets of First and Second Kind OPRL 

Our goal in this section is to prove (|1.13p through (|1.16p . We note that since 
decreasing iV to iV — 1 does not change the Poisson brackets for {ak}^Si, {bj}fSi^ 
nor the {Pj}fSi^, we can suppose n = iV in proving (|1.13p - (|1.16p but then increase 
N in the induction step. 

Proposition 2.1. ([TTSl) and ([LMjl holds for n = l. 

Proof. For N = 1, bi = c, Po{x) — 1, Pi{x) = x — c, and all Poisson brackets are 
zero. Thus, we need only show the RHS of (|1.14p is zero, that is. 



{x-c)-{y- c) 



1 







x-y 



which is obvious. 



□ 



Lemma 2.2. 



{alP„iw)}^-^alP,,^iiw) (2.1) 



Proof. Since 




Theorem 2.3. ([LT3)) and ([TU)) hold for all n. 
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Proof. Proposition 12.11 is the result for ri = 1. So, by induction, we can suppose 
(|1.13|) and ()1.14p for j < n and need only prove it for n + 1. By (12. 2p for j = n 
and the facts that {P„^i(x), P„_i(j/)} — {Pn{x), Pn{y)} = and P„ independent 
of On, Qn+i and P„_i independent of &„ and only the cross terms enter and 

-{Pn+i{x),P„+i{y)} {{x - 6„+i)P„(a;),a^jP„-i(?/)} - (a; ^ y) 

where (x ^ y) is like the first term but x and y are reversed. 
This is a sum of three terms: ^1,^2,^3 where 

2h = 2{x ~ bn+i)al{Pn{x),Pn-i{y)} -(x^y) 

'{Pn{x) - P„-i(y) - P„(y)P„-i(a;)) 



- y)an 



- Pn-l{x)Pn^l{y) 



(2.3) 



(2.4) 



(2.5) 



x-y 

by the induction hypothesis and the symmetry of {P„(y), P„_i(y)} under x y. 
Next, 

2t2 = 2{xPn{x),al}Pn-i{y) " ^ v) 
= al{x'' y)Pa-i{x)Pn-i{y) 

by Lemma [221 Finally, 

2^3 = -2{&„+i,a2}(P„(.T)P„_i(y) - {x ~ y)) 

= -al{Pn{x)Pn-l{y) - Pn{y)Pn-l{x)) 

There is a fourth term involving {6„+i, P„_i(?/)}, but this is zero. 

Clearly, ti^-t^+H =Q since (|2.5p cancels the first term in (|2.3p . and (|2.4p the 
second term. This proves that {Pn+i{x) , Pn+iiy)} — 0. 

Since {P„(a;), P„(?/)} = and {6„+i, P„(j/)} (since P„(y) only depends on 
{fljI^Cj), (El implies 

{P„+i(a;),P„(2/)} = {-a2p„_i(a;),P„(y)} 

= -a2{p„_i(a;),P„(y)} - P,,_^{x){al, P,,{y)} 

The first term is evaluated by induction and the second by Lemma 
The Lemma 12.21 term is +ia^P„_i(a;)P„_i(y) which cancels one term in 
{Pn-i{x),P^{y)}. Thus 



(2.6) 



2{P„+i(x),P„(2/)} = a2 



P„(y)P„_i(a;) -P„(a;)P„_i(y) 



y-x 



Pn{y){Pn+l{x) " {X - bn+i)Pn{x)) - {x ^ y) 



y-x 

iPn+liy)Pn{y) - Pn+l{y)Pnix)) 

X-y 



Pn{x)Pn[y) 



(2.7) 



proving the required formula for n + 1. In the above, ()2.7p comes from (12. 2p . □ 

The monic second kind polynomials, Qm can be defined as follows: We note that 
if J(6i, . . . , 6„; oi, . . . , a„_i) is the matrix (jl.ip with N = n, then 

P„(a;) = det(a; - J(6i, . . . ,6„;ai, . . . ,a„^i)) (2.8) 

as is well known. Q„ is then defined by removing the top row and left column, that 
is, 

Qnix) = det{x - J(62, . . . ,6„;a2, . . . ,a„_i)) (2.9) 
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Notice that J(6i, . . . , 6„; ai, . . . , a„-i) and J(6„, 6„_i, . . . , 61; a„-i, . . . , ai) are 
unitarily equivalent under the unitary USj = Sn+i-j, j — l,...,n. This shows 
making the a,b dependence exphcit (we include redundant variables; Qn is inde- 
pendent of 61 and ai): 

Pnix; 61, . . . , 5„; ai, . . . , a„_i) = F„(x; &„, . . . , 61; a„_i, . . . , ai) 

Qnix; bi,. . . , 6„; ai, . . . , a„_i) = P„_i(a;; 6„, . . . 61; a„_i, . . . , ai) 

Taking into account that reversing order changes the signs in (|1.14p and (|1.16p , we 
see that { , } changes sign, that is, (|1.14p implies (|1.16p . and we have 

Theorem 2.4. pTTS)) and pTT6)l hold for all n. 

Instead of using induction from P„,P„_i to P„+i,P„ and the above reversal 
argument, we can directly use induction from P„, Q„ to Pn+l^ Qn+i using 

Qn+i{x; &i, . . . ai, . . . , a„) = Pn{x; &2, ■ • ■ bn+i; 02, . . . ,a„) (2.10) 

and 

Pn+i{x; 61, ... , 6„+i; ai, . . . , a,i) 

= (a; - 6i)P„(x; 62, ... , 6n+i; 02, ...,»«)- a\Qn{x] &2, . . . , 6n+i; 02, ...,««) 

(2.11) 

that is, coefRcient stripping. 

3. Fundamental Poisson Brackets for OPRL 

Our goal in this section is to prove (|1.12p and, most importantly, (|l.lip . In the 
calculations below, we will compute Poisson brackets of functions of a, 6 and free 
parameter(s), x (e.g., P„(x)) and then set x to a value that is a function, h, of a, 6. 
The order of operations is important, that is, {., .} first, only then evaluation. We 
wiU denote this by {f,g}\^^h- 

Proposition 3.1. 

{P„(a;),P„(y)}|^^^^._y^^^ ^ {xj,Xk}W{xj - xi)'W{xk - x^) (3.1) 

In particular, 

(y)} = <^ {xj,Xk} = for all j, k (3.2) 
Proof. Since P is monic and its zeros are simple and (|2.8p holds, 

Pn{x)^l[{x-x,) (3.3) 

i 

Thus, by Leibnitz's rule, 

{P„(a;), P„(y)} = ^{xp,xj [](a; - x,) [](y - Xi) (3.4) 

p.q e^p l^q 

Setting X ~ Xj, y — Xk, all terms with p ^ j or q ^ k vanish and we obtain (|3.ip . 
([XT]) trivially implies in (IX^ and implies <^ in ((X^ . □ 

Proposition 3.2. FFe /laue t/iai 

n 

Qnix) ^^PjY\_i^ ^ ^k) (3.5) 
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Proof. The m-function is defined by 



pj 



3 = 1 ^ 



which is just 

TO„(z) = {5i, (J(6i, . . . ,fo„;ai, . . . ,a„_i) - zY^Si) 
which, by Cramer's rule and (|2.8p / (|2.9p , impHes 



m„(z) 



Thus, by 



(5n(2) = -m„(z)P„(z) 



which is (j3.5p . 
Proposition 3.3. 

{P„(a:),Q„(2/)} 



(3.6) 

(3.7) 
(3.8) 



□ 



{a;j,pfc} ]^(a;j - a;^) ]^(xfe - x^) (3.9) 
This implies that 

^ {^i,Pk} = i [SjkPj - PjPk] (3.10) 
Proof. Since we have proven that {xj, Xk\ = 0, p.3p and (|3.5|) imply that 

{P„(x), Q„(2/)} = - ^ {xp, pJ n(a; - a;^) X{{v " (^.ll) 

which immediately implies p.9p . 
In (fTTTC)) . we have, by (^3)1 . that 

Qn(a;)Qn(2/)L=^^,j,=:r, = {Pj/'fc}n(^J -a;£)]^(Xfc - ) (3.12) 

l^k 

The Bezoutian term in (|1.16p is if j ^ k, but is slightly subtle if j = fc because then 
the formal expression is 0/0. We thus take limits since, of course, {P„(x), Q„(?/)}, 
as a polynomial in x and y, is continuous. Since Pn{xk) — 0, for x ^ Xk, 

Pn{x)Qn{Xk) - Pn{Xk)Qn{x) 



X - Xk 



Wix-Xf) pk W{x-xt) 

i^k -I '-i^k 



Thus 



\im LHS of (IXT^ = Pk 

X^Xk 



Sjk 



(3.13) 



(3.14) 



Y]_{xk - xi) 

-i^k 

Thus (fTTC)) implies the right side of (|XTU)) . By (jXTTI) . the right side of (jXTUl) 
implies (frT6| . □ 

Remark. The symmetry of {a;j,pfc} under j ^ fc is equivalent to the symmetry of 
the right side of (|1.16p under x ^ y. 
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Proposition 3.4. Let j ^ k. Then 

_ PjPk PjPkPrn ST^ PjPkPjn 



Xk , . Xj Xfji Xk X 



Proof. f|3.15p comes from {Qn{x), Qn{y)} = 0. Explicitly, 

-{pj,Pfc}-RHs of jsm 



{Qn{x), Qn{y)}x=Xj, y=Xk 



Ue^jiX] - Xi) Ue^ki^k - Xi) 
by a straightforward but tedious calculation. (|3.16p then follows by using (|l.lip for 

{xp,pq}. a 

Remark. There is a sense in which implies (|1.12p since (|l.lip implies 

{Pn{x),Pn{y)} = which means {Qn{x), Qn{y) } = (since Qn is a F„_i for 
suitable Jacobi parameters), and this yields p.lSp . 

Finally, we want to translate the basic Poisson brackets relations (|1.14p / (jl.lSp 
and (jl.lip into assertions about the two-form which defines the symplectic struc- 
ture underlying the Poisson brackets. Recall a symplectic structure is defined by 
a two-form, w, that is, an antisymmetric functional on tangent vectors with the 
requirement that w is closed (which is equivalent to the Jacobi identity for the 
Poisson bracket; see, e.g., Deift [H]) and which is nondegenerate; that is, for all 
tangent vectors v ^ at p G M, there is at p so uj{v,v) ^ 0. Nondegeneracy 
implies d = dim(M) is even. 

Given such a two-form and function, /, the Hamiltonian vector field Hf is defined 

by 

u;{Hf,v) = {df,v) (3.17) 

for all tangent vectors v. Here df is the one-form '§i''^Xi as usual, and (■, ■) 

is the linear algebra pairing of vectors and their duals. 
The Poisson bracket associated to uj is then 

{f,g}^Hf{g) (3.18) 

Proposition 3.5. Let M be a symplectic manifold of dimension d — 2£ and suppose 
there is a local coordinate system {xj}j^i, {yjYj^i, so that 

e e 
u) = Wij dxi A dyj + Uij dxi A dxj (3.19) 

{i.e., no dyi A dyj terms) with U antisymmetric. Then 
(i) Nondegeneracy of uj is equivalent to invertibility ofthelxi matrix {Wij)i<i_j<i. 
We will denote its inverse by {W~^)ij. 

(ii) 

uj^ = i\ det(M^)da;i A dyi A rfa;2 A • • • A dye (3.20) 

(iii) 

{x,,Xj}^0 (3.21) 

(iv) 

{x,,y,} ^ {W~'),, (3.22) 
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Moreover, {yi, yj} — for all i,j if and only ifU = 0. 

Conversely, if a symplectic manifold has a Poisson bracket obeying p.2ip and 
([X^ . then uj has the form ((XTg)l . 

Remarks. 1. In (|3.20[) . uj^ means the ^-fold wedge product w A • • • A cj. 

2. W and U can be functions of x, y. 

3. These calculations do not use = (or, equivalently, the Jacobi identity). 
Proof. Write 



If t^ 0, pick s = W-^t and i = Y1=i ^k-^- If t = 0, r^ 0, and with s = 

we can take v = ^'^^af^' Thus, invertibility implies nondegeneracy. If W is 

not invertible, then det(VF) — and the calculation of (ii) shows uu is degenerate. 

(ii) is an easy calculation: using the distributive law to expand lo^ , any term 
with a dxi A dxj term has at least £ + 1 dx^s, and so is zero by antisymmetry. The 
only products of € dxi A dyj that are nonzero are of the form [(ix^(i) A A ■ • • A 
[(ia;^(^) A dy„(^^)\ with tt, ct permutations, and this is (—1)'''^ dxi A dyi A ■ • ■ A dy^. 
The sum over a for fixed tt yields det(VF) and then the sum over tt gives l\. 

(iii) , (iv). We begin by noting that if {Cj jjli is any coordinate system, if is a 
2£ X 2t antisymmetric real matrix, and 




21 




(3.23) 



is a symplectic form, then for 




(3.24) 



we have 




m— 1 



or (since ^rnk = -^km) 

au = [{-n)-^6.j\k 



(3.25) 



so that 



{QXk}^H^A^u) 
For our case if C = [xi, ...,X£,yi,.. .,yi), 



SI has the £ x £ block form 



-(SI ^)fcj 



(3.26) 




(3.27) 



for which 




(3.28) 



so (|33T|) and ([3:22| follow from 



10 



M. J. CANTERO AND B. SIMON 



For {yi,yj} = if and only if [/ = 0, we note {yi,yj} = if and only if 
W-^UW-^ = by (I3l8ll . 

For the converse, we note that p.2ip and p.22p plus antisymmetry imply fl^^ 
has the form (|3?28| for some U, and then |3?23| imphes ((3?T9)) . □ 

From this proposition, we can find the basic symplectic two-form in both (a, b) 
and {^J.,x) coordinates: 

Theorem 3.6. The symplectic form defined by ()1.3p / (|1.4|1 has the form: 

4 aJ^doiAdbk (3.29) 

i<e<k<N-i 

Remark. (|1.3p / (ll.Sp only define a two-form. That it is symplectic (i.e., duj = 0) 
follows from ((X29)) . 

Proof Since (O / ([L4l) say 

Proposition 13.51 fwhich. as we noted, does not use that the Poisson bracket obeys 
the Jacobi identity) says that (recall that fei, . . . , b^-i] ai, ■ ■ ■ , o-n-i are a set of 
coordinates) 

N-l 

LU — 

» J = l 

where {W~^Y has the form 



= da, A dbj (3.30) 



'ioi -iai ...\ 
ia2 -ia2 ... =i?(l-A^) (3.31) 

where D is the diagonal matrix j OkSki and iV=^o o i^is the standard rank 

iV - 2 nilpotent on M^-^. Thus 

W = [(£1(1 - N))-^Y = D-^{1 + N* + (iV*)2 + • • • + (iV*)"-i) (3.32) 

/4a-i ...\ 
_ 40^1 4a2"i ... , . 

~ 4ar^ 4a2-i 403^1 ... ^^'^^^ 

V / 

so ([5301 is (jX^ . □ 
Theorem 3.7. The symplectic form defined by (|1.3p / (11. 4p /las i/ie form: 

N 

2 ^ dxj A pj^ dpj + ^ fJy A dx^ (3.34) 
j=i ij 

/or some U . 

Remarks. 1. By X^jLi — c and X^jLi Pj = 1; the a;'s and p's are not independent, 
but they still define functions. 

2. One could compute U from (|1.12p . 
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Proof. Define for fc = 1, . . . , iV — 1, 

Vk = log (3.35) 
.Pn_ 

Then for j = 1, . . . , iV - 1, 

{xj,yk} = ^Sjk (3.36) 

since, by (jl.lip . 

{xj,yk} {Xj,\0gpk} - {Xj,\0gpn} 

= ih^jk-Pj) - i~Pj) 

{^j}f=i^: {yj}f=i^ are local (indeed, global) coordinates since X]j=i Pj — 1 ^^ts 
us invert p.35p . 

1 

so {xj,Xk} = and (|3.26p plus Proposition 13.51 imply 

JV-l 

u = 2 dxj A c?2/j + Uijdxi A (3.39) 
p.34p follows if we note that 

dVj = Pj^ dpj - pjj^ dpN 

and 



Pi = , ^jv-i „„, (3-37) 



= , ^jv-i 1 (3.38) 



Af-l 

dxj A dpN — —dxN A dpAT 



since X^jLi ^2;^ =0. □ 

4. Coefficient Stripping for Paraorthogonal Polynomials 

Paraorthogonal polynomials are defined [55] by 

Pn[z- {aJ}f^-o^/3) = z<i>Ar-i(z) - /3$^„i(z) (4.1) 

where /3 G i9D. Second kind paraorthogonal polynomials by 

QM{z-{a,}'^-^\p)=z^M-i{z) + ii-^*^_^{z) (4.2) 

where ^' are the second kind polynomials. These polynomials have been extensively 
studied recently [H [5l [181 Ull HI] and, in particular, the second kind polynomials 
are introduced in [111 HH] ■ The symbols P and Q are new but quite natural. These 
are relevant for the following reason: 

Proposition 4.1. Let 

N 

dp = ^p.jSz^ (4.3) 
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be a pure point probability measure on dD with each fj,j > 0. Let ^0,^1, ■ ■ ■ be 
the monic OPUC where, in particular, 

N 

■fN{z) = \{{z~z,) (4.4) 
Then the recursion relations define parameters {ctj^^S^ and aiq-i — f3 (z dD where 

N 

P={-ir+'l[-z, (4.5) 



Here 

Moreover, 



$w(z, d^l) = Pn{z; {«,};Lo, f3) (4.6) 



e'^ + z 



F{z,d[i)^ \ — d^i[z) 



Qiv(^;K}f^'o',/?) 

Piv(z;K}f."o':/?) 



(4.8) 



In particular, 



N 

In'-' 



Remarks. 1. Parts of this proof are close to results of Jones, Njastad, and Thron 
[22] as presented in Theorem 2.2.12 of [40]. 

2. (|4.9p is, of course, the analog of p.Sp . 

3. It is known that Imi^(e'^) is strictly monotone and pure imaginary between 
poles, so there is a single zero between poles. Thus, (|4.8I) shows the zeros of Pn and 
Qn interlace, a result of [HI US] obtained by other means. 

Proof. If d/i has the form (|4.3p . then L"^ [dD , dy) is A^-dimensional, so {^j}^"^^ 
span the whole space so 'i>Ar(e*^) — for a.e. 6. It follows that <I>Ar must obey 
(we use Pn for later purposes, even though we have not yet proven that $Ar is a 

Pn) 

N 

PN{z) = X{{z-Zj) (4.10) 



In particular, a^-i = — $Ar(0) is given by (|4.5p and so is a /3 G SB. Szego 
recursion thus says (|4.1|) holds. 

Using * for degree A^-polynomials, ()4.ip / (l4.2p imply 

P*N = -I3Pn Q*n = PQn (4.11) 

so 

Qn _ Qn 
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and (|4.8p is just (3.2.19) of 40J (this is proven for (ap, ■ • ■ , aw-i, 0, . . . ) with 
|aAr-i| < 1, but by taking aw-i ^ dD, one gets F = Q*n/Pn)- (SSI), and 

imply gH). □ 

Since we are about to consider changes in N, we shift notation from N to n. 
Since PmQn are defined via a boundary condition on q;„_i, if we want to prove 
something inductively, we need to "strip from the front," that is, remove ao. In 
[40l |4T| we called this coefficient stripping. For OPUC, this involves the formula 
for $„(z; {aj}"ro ) and <&* (z; ) in terms of ao and <i>„_i(z; {aj+i}Jrg^) and 

^n-ii^j {Q^j+i}j=o )■ Quite remarkably (and conveniently for use in Sections [5] and 
[6]), coefficient stripping for P„, Qn involves another P, Q pair even though P, Q are 
($,4") mixed! It turns out simpler to state things in terms of 

C„(z;{a,}^^o^/3) = ^^i±^ (4.13) 

^n(^;K};Co^/3) = ^^^^ (4.14) 

which we name analogously to cosh and sinh. 

Theorem 4.2. Let f3 e dD be fixed and {a^}"^^^ e B"-^. Let C„ denote 
Cn{z]{ajYjZQ,l3), and similarly for Sn- Define Cn~i by 

Cn-i = Cn-i{z; {a.+ij'JllP) (4.15) 
(with ao removed), and similarly for Sn-i- Then 

Cniz) = ziCn-iiz) - aozSn-i{z)) (4.16) 
Sn{z) = -aoCn-i{z) + Sn-liz) (4.17) 

(4.18) 



Proof. As usual, let 



A(a,z)=( 1 " (4.19) 



—az 1 

so with (5+ = (p) e C^, S- = d*) e C^, Szego recursion says 



n— 1 



n — 1 



A(a„-2, ^) . . . A(ao, ^)(<^+ + S-) (4.20) 
A(a„_2, ^) . . . A(ao, ^)(<^+ - <5-) (4.21) 



Let B{p) = ( 

0-/3) ^^"4 conclude 

P„ = ((5++(5_),(B(/3)A„_2(a„_2,z)...A(ao,;2))(<5+ + <^-)) (4.22) 

Q„ is similar with the rightmost ((5+ + (5_) replaced by ((5+ — (5_). Thus, C„,S'„ 
correspond to using (5+ and (5_ , and we have that 

(^^^ - A(ao, z)*A(ai, z)* . . . A(a„_i, z)*i?(/3)* (4.23) 

which implies 



Aiao,zr[ (4.24) 
which is ((4A6)) -||4T7 |1 . □ 
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Remark. One might think the transposes in (I4.23P should be adjoints, but the shift 
of d± from the right to the left side of the inner product introduces a complex 
conjugate. 

While we won't need it, we note the induced formula for P„, Qn- 

Theorem 4.3. //P„(z) = P„(z; {a,};^o^ /3) and P„_i(z) = Pn{z; {a,+i}]I^ , P), 
and similarly for Qn, Qn-i, then 

Pn^ ^{z -ao + l- aoz)Pn-i + i (z - ao - 1 + aoz)Qn^i (4-25) 

Qn = 5 (z + ao + 1 + aoz)Qn~i + 5 (z + ao - 1 - aoz)Pn-i (4.26) 

Proof. This is immediate from (|4.24p . together with 

,QnJ V 1 -1 / \Sn 
'Cn.i\ ^1(1 1\ (Pn-l 



and a calculation of 



1 1 \ X / 1 1 

1 



^ lA'(ao,z) ^ _^ 1 □ 



The algebraic simplicity of the recursion of (C„, Sn) relative to (P„, Qn) is rem- 
iniscent of Schur vs. Caratheodory functions (see Section 3.4 of '^). This is not a 
coincidence. Note that by 



Cn ^ Pn + Qn ^ 1 - [{-Qn)/Pn] ,^ . 

Sn Pn -Qn 1 + [{-Qn)lPn] ^ ' ' 

1 - F 



l + F 



m = (4.28) 

Note that Sn is a polynomial of degree n — 1 and Cn is a polynomial of degree 
n vanishing at zero, so (|4.28p is a ratio of polynomials of degree n — 1. —Sn and 
z^^Cn are thus essentially para-versions of the Wall polynomials. Indeed, by the 
Pinter-Nevai formula [40l Theorem 3.2.10], 

z-'CJz) = zB:_2{z) + (3An^2{z) 

Sn{z) ^ -ZAI_^{Z) - pBn-2{z) 

5. Poisson Brackets of First and Second Kind OPUC 

Our goal in this section is to prove (jl.lSp and (|1.19p . Once we translate to C, S, 
the induction will be even simpler than in the OPRL because C„_i, Sn~i are only 
dependent on {aj+i}"Jo ^'^^ tiave zero Poisson bracket with ao, so there are no 
analogs of terms like (j2.ip . Since we are varying the index, we use n in place of N . 
We will begin by showing p.lSp / p.l9p is equivalent to 

{Cn{z),Cn{w)} = = {Sn{z),Sn{w)} (5.1) 
sn ( \ Q ( \\ ■( Cn{z)wSn{w) - Cn{w)zSn{z) \ , 

{6„(z),5„(w)} = -I (5.2) 
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Proposition 5.1. ()1.18p and the symmetry o/{P„(z), Qniw)} under z <-> w implies 
(f5?T|) and 

{C„(z), S„{w)} = -i {P„{z), Q„(w;)} (5.3) 

Similarly, (|5.ip and the symmetry of {Cn{z), Sn{w)} under z ^ w implies (|1.18|) 
and (lO). 



Proof. Given the antisymmetry of the Poisson bracket, the symmetry of 
{Pn[z),Qn{w)} can be written 

{Pn{z), Qn{w)} + {Qn{z), P„H} - 

from which (|1.18p implies (|5.ip . and conversely, (|5.ip and the symmetry of 
{Cn{z),Sn{w)] implies ^M- 

The antisymmetry plus (jl.lSp also shows 



{c„(z),5„H} = i2{g„(z),p„H} 

which is (15.31). and the converse is similar. 



□ 



Next, we do a calculation showing the equality of the left side of (|5.2p and of 
(fTlQl) using P„ = C„ + Sn and Q„ = C„ - S'ni 



RHS of (fLlQl) 



[(r 


+ w\ 




— w J 



z + w 



'-2Cn{z)Sn{w)+2C„{w)Sn{z)] 

+ 2C„(z)5„(u;) + 2C„H5„(z) 

z + w 



1 ^C„(z)5„(w 



2Cniz)wSn{w) - 2Cn{w)zSn{z 



+ l'>Cn{w)Sniz) 



z — w 
= -2[RHS of (lO) ] 



(5.4) 

consistent with (|5.ip . 

Thus, (jl.lSp / p.l9p is equivalent to (|5.ip / (|5.2p . which we proceed to prove in- 
ductively, starting with n = 1. We let G„(z, w) denote the right side of (|5.2p . 



Proposition 5.2. We /laue t/iai 

Ci(z) = z 5i(z) = -^ 
T/iws, 01(2,11;) = and (|5TT|l . (|5?2|l /lo/d for n = l. 
Proof. Since 

Pi(z) = z-^ Qi(z) = z + ^ 
(15.51) is immediate, and thus. 



(5.5) 



(5.6) 



Gi(z, w) = i/3 



zm; — wz 

z — w 



= 



Since Ci, iSi are independent of {ctj}, all Poisson brackets are zero and (|5.ip . ([5 
hold for n = 1. 



□ 
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Remark. It is an interesting exercise to prove 

z vu 

(Pi (z)Qi M - Pi iw)Qi (z)) Qi{z)Qi{w)+ Pi (z)Pi {w) 

z — w 

is from (j5.6p . and explain why —QQ appears with a minus sign. 
Lemma 5.3. G„ and G„_i are related by 

Gn{z, W) = ZplGn-l{z, w) - iplzSn-l{z)Cn-l{w) (5.7) 

Proof. By (|iT^ . (|iT7)) . and the fact that 

m/,g)(^,^)= ^^^^-'^^"^~^^"^"^^^ (5.8) 
z — w 

is zero if / = g: 

G'„(z,w) = -iB{Cn,XSn){z,w) 

= —iB{XCn-i — a^XSn-l-, -agXCn-l + XSn-l) 

= -z[S(XC„-i,^5„_i) + |aopS(X5„_i,XC„_i)] 

= -iplB{XCn^l,XSn-l) (5.9) 

where we use X f for the function 

{Xf){z) = zf{z) 

and (|5.9p using antisymmetry of B in / and g and 1 — |ao|^ = Po- 
Next, note 

zf{z)g{w) - wf{w)g{z) 



B{Xf,g){z,w) 



z — w 

f{z)g{w) ~ f{w)g{z)\ (z-w) 



f{w)g{z) 



z — w 

= zB{.f,g){z,w)+ g{z)f{w) 

to conclude 

- iB{XCn-i,XSn-i) = zG„_i(z,w) - izSn-i{z)Cn-i{w) (5.10) 

([5:91 and ((5J0)) imply dSJ]). □ 

Theorem 5.4. jEI]) -([52]) hold for all n. (fTTSl) - (frTg]) hold for all n. 

Proof As noted after Proposition[5lI]and (fO]) . (I1.18|) . (I1.19|) is equivalent to ([5?T|) 
and (j5.2p . So, by Proposition 15.21 and induction, we need only prove 

(lOjl - ljO)) for n - 1 dSHl-jlSl) for n (5.11) 

By ()4.16p and (I4.17p . {q!o,q;o} — {ao,ao} = the symmetry of G„_i(z,w) in 
z ^ w, and the independence of S'n-i, G„_i on ao, aoj we have 

dm])- ([521) for n - 1 ^ jSU) for n 

(Notice how much simpler this is than the analog for OPRL!) 
On the other hand, by (lilB and (|4l7l) . 

{C„{z),S„{w)} = z{G„_i(z) - aoSn-iiz), -aoG„_i(w) + S'„_i(w)} 

= zplGn-i{z, w) + {ao, ao}zSn-i{z)Cn-iiw) 

= zplGn-i{z, w) - ipo2'S'„_i(z)G„_i(w) (5.12) 
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= Gn{z,w) (5.13) 

([5?T2ll uses (fTTTl) and ([5?T3l) follows from (fST]) . We have thus proven ([O]) for 
n. □ 



6. Fundamental Poisson Brackets for OPUC 

In this section, we will establish (|1.23p . Let {zj}^^^ be the zeros of Pn{z), Oj 
given by (|1.2ip and /ij by (|1.20p . Unlike the OPRL case where xi < ■ ■ ■ < x„ 
allowed global variables, these are only defined locally, both because of the z's not 
having a unique initial zi and because of 27r ambiguities in 9j. 

Throughout this section, we fix /3 e 9D and consider the 2n — 2-dimensional 
manifold D"~^ of Verblunsky coefficients {aj}^^^ and a„_i = f3. The associated 
finite CMV matrix C has 

det(C) = (-1)""V 
The associated spectral measure given by (|1.2ip thus has 



ri-l 



and, of course, 



^0,=arg((-i; 

n 

so as local coordinates {dj}"Zi U {fJ-j}"'Zi coordinatizes 



■f.n-1 



(6.1) 
(6.2) 

(6.3) 
(6.4) 

/? never appears 



explicitly below — this is not surprising since a rotation can move /3 to 1. 

As usual, one can replace ReQ;,ImQ; by "independent" coordinates a, a with 



da = d(Re a) + i d(lni a) da — d(Re a) — i d(Im a) 



d_ 

da 



d 



1 d 



d_ 

da 



d 



1 



d 



d Re a id Im a 



d Re a id Im a 

Thus, the Poisson bracket defined by (|1.17p can be defined in general by 
U l-Y^- 2/^ df dg df dg 



daj duj 



We will denote 



(f'aj — d{Reaj) A d{lmaj) 



— daj A daj 



Proposition 6.1. For OPUC, 

{Pn{z),Pn{w)}\ 



{9k,e,}Y[izk-z,)Y[{2 



In particular, 



{Ok,Oj} = 



(6.5) 
(6.6) 

(6.7) 

(6.8) 
(6.9) 

(6.10) 

(6.11) 
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Proof. Since Pn{z) — YG=ii^ ~ -^^)' have 

LHS of (Eini) = {zk, zj} Y\{zk - zi) - zt) 

(|6.10p follows if we note 

(fTTSl then iniphes (|OT|) . □ 
Proposition 6.2. For OPUC, 

{Pn{z),Qn{w)}\^^^^^^^^^ = ~i{9j,lXk}2ZjZkW_{z, -zi,)W{zk -zi) (6.12) 

while 

RHSof mBL=z, 

(6.13) 



-{2zki-ikSjkizj + Zk) - 2zj2zkfijtJ.k}Y\_i^j ~ ^i) II(^fe 



e^j i^k 

In particular, 

{6'j, ^ik} = ^^]6jk - Aij/^fc (6.14) 
Proof. Given (|4.9p . setting z = Zj, w — Zk picks out the {zj,fik} term, that is, 



LHS of (l^ia = {-Zj,Aife} 



]J(zj-z<,) 2zkY[{zk-ze) 

e^j i^k 



where 2zfe is (w + Zfc)|^^^^. Since 

{-Zj,fik} = -izj{0j,fik} 

we obtain ()6.12p . 

Setting w = Zfc in the right side of p.l9p gives (using P{zk) — 0, (3(z/c) 

/ife(2zfe)n^_^j^(zfe - 0^)) 

RHS of ^ML^z, = -'t^I^^Mz + Zk) - Q{z)fik](2zk) l[{zk - z,) 

) i^k 



Since P{z)/{z ~ Zk) — Yie^ki-^ ~ ■^^)' &rst term is zero as z Zj if j ^ fc, and 
otherwise is JJe^ji^j ~ ^-d- This yields (pl^ . (pT^ plus ((O^ yields (|5Tl|l . □ 

From {(5(z), Q(it;)} = 0, one can derive the formula for {/ij,/!^}. We do not 
provide details. 

Finally, for this section, we compute the symplectic form, w, in da, and dQj, dfij 
coordinates. 

The da coordinate is especially simple since in those coordinates, the Poisson 
bracket is a product structure. By p.l7p and Keaj — \{aj +Q;j), Imofj = ^{uj — 
aj), we get 

{Reaj,Imafc} = — 
so by the U — Q case of Proposition 13. 51 (and (|6.9p ): 
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Proposition 6.3. The symplectic form for OPUC with a„_i — 13 ^ 9D fixed is 

(6.15) 



= — d(Re aj ) A d(Im aj ) 

n-2 . 



-0 Pi 



(6.16) 



The calculation for 9, jj, coordinates is essentially the same as for OPRL with x, p 
coordinates since '}2^=i is constant. Since the ^ in (jl.lip is missing from (|1.23p . 
the 2 in (|3.34p is absent and we get 

Proposition 6.4. The symplectic form for OPUC in terms of {fij,9j} is 

N 

" 9k (6.17) 



w = ^ ddj A fj.j ^ dfi-j + ^ Ujk ddj A ddk 



for some U. 



7. Application 1: Jacobians of Coordinate Changes 

This is the first of three sections in which we present applications of the fun- 
damental Poisson brackets prTTj) for OPRL and (fOS]) for OPUC. Those in this 
section and the next are not new, but we include them here because the literature 
on these issues is not always so clear nor is it always emphasized that all that is 
needed are the fundamental Poisson brackets. 



Given two local coordinates xi, 
defined by 



and j/i, . . . , their Jacobian 



\dy\ 



named because 



d{x) 
d{y) 




det^ 


dxA 
dvk) 


d''x 




d{x) 
d{y) 





(7.1) 



(7.2) 



relates the local volume elements. Of course, differential forms are ideal for this. 



dxi 



Yj ^ dyk means 



dxi A • ■ • A dxr^ 



detfe 

\oyk 



dyi A • • • A dyk 



(7.3) 



Our goal here is to compute the two Jacobians | | for OPRL and | ^^^g(g'^™"^ I 
for OPUC. It appears that these were first found for OPRL by Dumitriu-Edelman 
[TT] . and motivated by that, for OPUC by Killip-Nenciu [28]. Their argument 
was indirect and Killip-Nenciu asked if there weren't a direct calculation. This was 
provided by Forrester-Rains |15j using forms and the continued fraction expansions 
of resolvents. Their argument is not unrelated to the one below, but is lacking the 
connection to Poisson brackets and we feel is more involved. In ^29J, Killip-Nenciu 
remark that Dcift (unpublished) explained to them how to go from the symplectic 
form for Jacobi parameters rewritten in terms of {x,p) to | g|°|p| | and they then 
do the analogous OPUC calculation. As they remark, this is buried at the end 
of a long paper, and the fact that only the form of Poisson brackets is involved is 
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obscure. We present this idea here to make the OPRL argument exphcit and to 
present the OPUC argument without a need for Lie algebra actions. 

Theorem 7.1. On the 2N — 2 -dimensional manifold of Jacobi parameters with 
X)jLi fi'Xed, use {o.j,bj}^Si^ and {xj, Pj}^Si^ as coordinates. Then 



d{a,b) 




■ 2 


d{x,p) 







nN 



(7.4) 



// one considers the 2N — 1-dimensional manifold without X^jLi fixed {equiv- 
alently, without X^j^i fixed), then (|7.4[) is still true, but where the Jacobian 
is from (ai, . . . , a^-i, 6i, . . . , bj^j) to (xi, . . . , xn, Pi, ■ ■ ■ , Pn-i) rather than from 
(ai, . . . ,aN^i,bi, . . . ,bN-i) to {xj, . . . ,XN-i,Pj, ■ ■ ■ ,Pn-i)- 

Proof. Let lu be the symmetric form on the 2Af — 2-dimensional manifold determined 
by ^M-^M- By (15^29]) . the {N - l)-fold wedge product 

w A • • • A w = 4^"^(7V - 1)! Yl aj^ dai A dbi A ■ ■ ■ A daw-i A dbN-i (7.5) 

By p.34p , noting all terms with at least one dxk A dxj must have a repeated x 
and so vanish, 



to A - ■■ Alo 



N 



2^-\N -IV-Y. f\ Pk\dxk A dpk) 

r N 



dxi A dpi A ■ ■ ■ A dxN-i A dpN-i 



N 



k=l 



2^-i(7V - 1)! Jl pl^ \dxi A dpi A • ■ ■ AdxN-i A dpN-i 



(7.6) 
(7.7) 



Here (|7.6|) is obtained from 
dxN - 



(7. 



- dxj dpN = — dpj 

and (|7.7p uses 

N N N N 

T.UPk'-UPk'T.p^-UPk' 

j=i k=i j=i k=i 

l|7.4p is immediate from the equalities of the right side of (|7.5p and (|7.7p . 
To get the results for the 2N — 1-dimensional manifold, we note that for each 
fixed ^j: have 

RHS of (1731) = RHS of (I73D (7.9) 

Since X^jLi = SjLi > can wedge the right side of (|7.5p with dbi + • ■ • + d^Ar 
and of (|7.7p by dxi + • ■ • + dx^v and so obtain the analog of (|7.9p but with the extra 
variable. □ 
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Theorem 7.2. On the 2N — 2 -dimensional manifold where un-i = (i ^ 91D) is 
fixed and 

Uj ^ Uj + i Vj 

(w, u e M), use {uj,Vj}^SQ^ and {^jiMjl^L^^ '^^ coordinates. Then 



d{u, v) 



(7.10) 



If we consider the 2N ~2- dimensional manifold without fixing f3 and use coordinates 
{uj,Vj}fSf^ U {V'} where (3 = e*'^ and {%}{Li U {iJ-j}^^i , then 



d{u,v,Tp) 



Proof. By ((05| . the iV - 1-fold product 

• A-2 



RHs of (rrruii 



(7.11) 



w A • • • A 



= {N- 1)! 2^-1 n pT^^ du^ A d«o 



A • • • A (i?iAr_2 A (iuAr_2 (7.12) 



where, by (|6.17p and the same calculation that led to (|7.7p (i.e., (|7.8p with p replaced 
by A*), 



JV 

a; A • • • A ti; = (A^ - 1)! rf^i A d^i A • • • A dO^-i A d^^-i 



(7.13) 



l|7.10p is immediate and (|7.1ip then follows from 



N 



di: = Y^ de.j 
i-1 



on account of 



□ 



Note that in comparing (|7.4p and (|7.10p with Forrester-Rains [TF and Dumitriu- 
Edelman jTT], one needs to bear in mind that their qj are related to pj (and ^j) 

by 



and thus 



and therefore 



Pj = 1j 



dpj 2dqj 



Pj 1] 



dpi... dpM-i _ ^Tv-i J_ dqi... dq^-i 



Pi ■ ■ - Pn 



Qn qi ■■ -Qn 



(7.14) 
(7.15) 

(7.16) 



so they have no 

2-(^-i) in dZH) but have a 2(^-i) in HjTJU^i . There is an extra 
factor of l/gAT (which is in [TS] but missing in [TT] due to the fact that their d^~^q 
of [TT] is the measure on a sphere, not EucHdean measure). 
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8. Application 2: Generalized Toda and Schur Flows 

In this section, we want to show how (jl.lip aUows the exphcit solution of flows 
hke (|1.6p and (|1.7p in the p (resp. p) variables. As mentioned earlier, these results 
are not new. Theorem 18.11 appears in [9] and many times earlier. Theorem 18.31 
appears already in [291 as their Corollary 6.5. 

Theorem 8.1. Let f he a real-valued function on R and let 

N 

ff-E/(^^) (8-1) 

and let a, b solve (with ' — d/dt) 

dk^{H,ak} bk^{H,bk} (8.2) 

for some initial conditions for the Jacobi parameters. Then {o,k{t),bk{t)}^^^ are 
the Jacobi parameters associated to 

X, it) = X, (0) p, it) = Pj\ \ (8.3) 

Remarks. 1. Since all that matters is / at the xj, we can restrict to polynomials 
where if f{x) = a;™, then H = Tr(J'"). 

2. {xj,Xk} — implies all Tr(J™) are constants of the motion. 

Example 8.2. The Toda flow (HI]) / ILT]) corresponds to H = 2Tr(j2), so f{x) = 
2x'^ and f'{x) = Ax, and (|8.3|) becomes 

which is well known to solve Toda [3H [S] . 

Proof. The flow is generated by the vector field {H, ■ } so 

x,^{H,Xj} p,^{H,p,) (8.5) 

Since {xj, Xk\ — 0, (|8.5p implies Xj = 0, that is, Xj{t) = Xj{{)). 
By (inil),ifj^A^, 

Xk,— \^ PN^ixk^Pj} - -3- {xj^Pn} 
PN J Pn 

= 5 PN^^kjPj - PjPk)- ^ PN^iSkNPj ~ PjPk) 

= \ (4, - ^kN) (8.6) 
Pn 

Thus 

^l^\-\\r{^i)-!\xN)\^ 
dt\pN J Pn 

Defining yk{t) by p.35p . we get 

yk{t)=yk(S)) + lt{f\x,)^ f'{xN)) 
which by ([3:371 yields ((831). Q 



7) 
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Theorem 8.3. Let f be a real-valued function on 9D 

AT 

F = ^/(e^«.) (8.8) 

and let a solve 

<^j^{H,a,} (8.9) 

with some initial conditions |aj(0)| < 1 and boundary conditions a_i = — 1, 
ajv-i = (3. Then {aj{t)}^Si^ U {aN-i — P} are the Verblunsky parameters as- 
sociated to the N -point measure with parameters 



-it) = 6,(0) ^^,{t) = ^^^^ > (8.10) 



°*g(^"'^ Vj(o) 
Etie*^(^"'=Vfe(o) 

where g is the function given by 



9ie'') = -^f{e^') (8.11) 
Example 8.4. If fie'") = 2 sin 6*, then H = \Tt{C - C~^) and is the Schur 

flow [D 1 m uni 1^ 

= - aj-i) (8.12) 

e" 



g is 2 cos 6 and (|8.10p becomes 

„2tcos(e,)^^.(0) 



M.W- ^^ (8-13) 



Proof. As above, = implies 6j{t) = 9j{0). The calculation oi ^j{t)/ ^^{t) 

is identical to the one in Theorem 18.11 except there is no factor of i and df{x)/dx 
is replaced by d/(e*«)/d6'. □ 

By taking suitable limits, we get that certain measures on R (with all moments 
finite) as initial conditions for difference equations on {aj,bj}°^i or {aj j^i can 
be solved by 

Jetf'(-) dpt^x) 
See, for example, the discussion in [42] and references therein. 



9. Differential Equations for OPRL 

We have just seen that symplectic flows are naturally deflned on Jacobi parame- 
ters. On the level of measures, they are solved by (|8.3p . Here, we want to consider 
the differential equations on OPs induced by these flows. In terms of (|8.3|1 . we will 
suppose / is a polynomial; explicitly, 

k 

y'{x)^J2^,x^ (9.1) 

We drop the cq term since it drops out of dpt- Here is our main result: 
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Theorem 9.1. Under the map, 

with f given by (j9.ip . the monic orthogonal polynomials Pn{x\ dpt) — Pn{x] t) obey 

k k / n \ 

n+l,n — ^+1 

n aAPn-i{x;t) (9.3) 

£=1 j = l ^j=n-e+l ^ 

where J{t) is the Jacobi matrix ofdpt- 

Remark. The vector indices associated to J start at 1, so (with pj — Pj/\\Pj\\) 

Jjk = (P]-i,xpk-i) (9.4) 

■^n ^ ^] ^ (Pj-i^xp-j-i) (9.5) 

Jjj+i = Jj+i 3 = aj {pj,xpj^i) (9.6) 

Example 9.2. Take ■^f'{x) — x (the Toda case). Then J„+i,„ is aU that is relevant, 
and it is a„, so (|9.3p becomes 

Pn = -alPn-i (9.7) 
This Toda case is known; see, for example, Peherstorfer ^38J and Barrios-Hernandez 

m- □ 

Example 9.3. Take |/'(x) ^ x^. We have 

(>^^)n+i,Ti-i = a„a„_i (J^)„+i^„ = a„(6„+i + b„) 

so 

Pn = -al{b„+i + b„)Pn-i - alal_iPn-2 (9.8) 

□ 

Remark. Using (P„, P„) — a„a„_i . . . ai and (|9.5p . one can easily go from differen- 
tial equations for P„ back to those for a„ and b„. 



Proof of Theorem \9.1\ Since P„ is monic, P„ is a polynomial of degree n — 1. Thus 

n 

Pn = ^{Pn,Pn-i)Pn-e 
1=1 

n 

= -Y,{Pn,hf'Pn-l)Pn-l (9.9) 



RHS of (El 



where (j9.9p comes from 

O = — / PnPn^e d^t 

= {Pn,Pn-e) + {Pn,Pn-e) + {Pn, ^ f'Pn-t) 
= {Pn,Pn-t) + {Pn, \ f'Pn-l) 

since Pn-i is a polynomial of degree n — £, and so orthogonal to P„. □ 
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One can also ask about derivatives of pn = Pn/\\Pn\\- Since 



Pn 



Pn 

\Pn\ 



Pr. 



\\Pn 



\P, 



\Pn\ 



(9.10) 



we clearly need only find ^ log ||P,i||. Following [HI Sec. 9.10], 
Proposition 9.4. Under the map, 

with f given by (j9.ip . the monic orthogonal polynomials P„(a;; dpt) — Pn{x] t) obey 

k 



I log||P„|| ^ i5]c,[(J^)„+i „+i - Unii] 



Proof. Let 



dpt 



Since Pn is orthogonal to P„, 



since 



Pn{h.ndpt = \\Pnr{Pn,anpn) 



and we have (|9.4p . Similarly, 



= (Pl^ for n = 



Since 



|log[l|P„fiV*]-|logiV. 



(9.11) 

(9.12) 

(9.13) 
(9.14) 

(9.15) 



(Pl^ and (Pl^ imply (PTTjl . 

By (p?3)) . ((9101) . and ((9?TT|) . wc immediately have 
Theorem 9.5. Under the map, 

wzi/i /' given by (|9.ip . ^/ie orthonormal polynomials pn — Pn/\\Pn\\ obey 

, k -K k r- k 

Pn = -Pn[^'^Cj[{J^)n+l n+1 " {J')l l] ) " X] X] (^^^ 



(9.16) 



□ 



«=1 ^-3 = 1 



Pn-l 
(9.17) 
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10. Differential Equations for OPUC 

For OPUC, one looks at flows generated by Laurent polynomials, f{z), real on 
dB. If 

j=-e 

where &_j — bj, then 

The analog of the P„'s are the unnormalized CMV and alternate CMV bases 

Yn ^ PO ■ ■ ■ Pn~lXn (10.3) 
Xn = PO ■ ■ ■ Pn-lXn (10-4) 

where {Xn}^=o and {a;„}5^g are the CMV and alternate CMV bases (see [101 
Sec. 4.2]). Then the exact same argument as led to Theorem 19.11 leads to 

Theorem 10.1. Let dpt be given by (|10.2p with g given by (|10.ip . Then ifYn{t) = 
Yn{dpt), 

m— k——i ^ j—m ^ 

with the same equation for Xn except the (C'^')m„ is replaced by {C'^)nm- 

Remark. These are analogs of Proposition [HH] and Theorem 19.51 

Example 10.2 (Schur Flow). We have 

g{z) = z + z-^ (10.6) 

Ismail f3P (see also Golinskii [TO]) obtained the following differential equation for 
the monic OPUC: 

$„ $„+i - (z + a„a„„i)$„ - pl_j^^n-i (10.7) 

We want to show this is equivalent to (|10.5[) . From our point of view, (|10.7p is 
unusual since our arguments show that $„ is a linear continuation of {^j}"ZQ 
while (|10.7p has <i>„+i and <&„. (Since there is a z$„, in fact, one can hope — and it 
happens— that RHS of pUT)) is Y.''j=q Cj<^j-) 
Since $„+i = z$„ - an^*„, 

RHS of (UnZl) = -a„{K + an-M - p'_i*n-i (10-8) 

Since (see (1.5.41) of gOj) K + - pl-iK-i, 

RHS of (Unjl) - -ap2_i$;_i - (10.9) 

verifying that RHS of (|10.7p is indeed a polynomial of degree at most n — 1. 

(|10.9p mixes Vs and X's, so we use (see (1.5.40) of piT) $„_i = pl_2Z^n-2 - 
an-2*;!;_i to find 

RHS of (Unil) = -Pl-lian - an.2Wn-l - Pl-lPl-2{-^^n-2) 

which one can see is exactly (|10.7p or its X„ analog (depending on whether n in 
$„ is even or odd). □ 
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11. Application 3: Poisson Commutation for Periodic Jacobi and 

CMV Matrices 



In this section, we want to show how the fundamental relations (|l.lll) / (|1.23p on 
finite Jacobi and CMV matrices yield a proof of the basic Poisson bracket relations 
for the periodic case. For OPRL, the original proofs of this relation by Flaschka [14] 
via direct calculation of Poisson brackets of eigenvalues is simple, but for OPUC, the 
two existing proofs by Nenciu-Simon and Nenciu [34|, I35j are computationally 
involved. 

Given , G or {ajl^Zp G C (in the OPUC case, we assume p is even 

for reasons discussed in Chapter 11 of |4lj), we define periodic Jacobi and CMV 
matrices by first extending the parameters periodically to Z, that is 

Gj+p^Qj b-i+p^bj aj+p=aj (11.1) 

and then letting J, C act on (Z) as two-sided matrices with these periodic param- 
eters (two-sided CMV matrices are discussed on pp. 589-590 of [41]). 

In each case, J, C act on €°°(Z) and commute with : {u„} {w„+p}. For 
each e'^ £ dU, J, C leave invariant the p-dimensional space 

l'^ = {u I S^u = e^'^u} 

and so define J{6),C{0) as operators on i"^ . In terms of the natural basis {S^}^^i, 



J{9) has a matrix of the form 

( 

ai 



k ^ j mod p 



(11.2) 



ai 
62 02 



\e-'^ap 



e''ap\ 




bp J 



(11.3) 



of a Jacobi matrix with an extra term in the corners. The C{0) are described on 
pp. 719-720 of lilj. 

In each case, there is a natural transfer or monodromy matrix, Tp{z) on C^, so 
if w G C°° solves (J — z)u = (resp. (C — z)v = 0), then 



Constancy of Wronskians implies 

det{Tp{z)) 



Up+l 
Ur, 



(OPRL) 
(OPUC) 



(11.4) 



(11.5) 



(11.6) 



One defines the discriminant, A(z), by 

^ \Tr{Tp{z)) (OPRL) 
\z-P'''Ti{Tp{z)) (OPUC) 

(recall p is even in the OPUC case), z is clearly an eigenvalue of J{9) (resp. C{6)) 
if and only if e*^ is an eigenvalue of Tp{z) (resp. z"P/^Tp(z); see |4T1 p. 722]). So, 
by (mU), 

z is an eigenvalue of J{e) or €{9) ^ A{z) = 2cos(6') (11.7) 
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For later purposes, we want to note that in neither case is ^{z) automatically 
monic. Rather, 

f zP + lower order (OPRL) 

XjTJ; (^''^^ + ^"^^^ + orAexfi in between) (OPUC) ' 

(see [45] and Chapter 11 of gl]). 

W X (0, oo)''"^ (resp. W) are given Poisson brackets by 

{6fc,afc} = -iafc /c = l,2,...,p (11.9) 
{6fc,afc_i} = iofc-i k = 2,...,p (11.10) 
{6i,ap} = iap (11-11) 

and 

{a-,, Ckfe} = {Q;j,aA;} = {ckj, afe} = -i p| Jj7c j, /c = 0, . . . ,p - 1 

One big difference between OPRL and OPUC is that in the OPUC case, the 
center of the Poisson bracket (i.e., the functions that Poisson commute with all 
functions) is trivial, that is, {■, ■} defines a symplectic form. But in the OPRL case, 

p p 
Y,b, and (11-12) 

are easily seen to have zero Poisson brackets with all a's and &'s, so the two-form 
defined by {■, ■} is degenerate, but on the subspace where 

p p 

bj = P and JJ^ aj = a (11.13) 

the form is nondegenerate and the Poisson bracket defines a symplectic form. 
Our main results in this section are: 

Theorem 11.1 (Flaschka [H]). For OPRL: 

{A(x),A(y)} = (11.14) 

If \j(P) are the simple eigenvalues of J (9), then 

{X,{9),Xk{9')}^0 (11.15) 

Remarks. 1. (|11.15p is normally only stated for 9 = 9', but we will see it is 
immediate from (|11.14p and (|11.7p . 

2. For 7^ 0, TT, all eigenvalues are simple and Xj{9) are global nonsingular func- 
tions. For = 0, TT, there are subvarieties of codimension 3 where some eigenvalues 
are degenerate. Xj{9) are smooth away from this subvariety and (jll.lSp only holds 
there. 

Theorem 11.2 (Nenciu-Simon [37]). For OPUC: 

{A(z), A(w)} = (11.16) 
If Xj{9) are the simple eigenvalues ofC{9), then 

{X,{9),Xk{9')} = Q (11.17) 
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and 

\X{pk,Xm\^Q (11.18) 

fe=0 ^ 

Remarks. 1. As in the OPRL case, \j{9) are simple \i 6 ^ Q^ir and of! a closed 
subvariety of codimension 3 if = or = tt. 

2. 0^=0 inverse leading term in A(^) for OPUC. The analog for OPRL 

is Hfe^i ^k- Of course, {0^=1 ^fc, = 0, but we don't say it exphcitly since 

for OPRL, Jlfe=i '^fe Poisson commutes with any function! 

We note two general related aspects of these theorems. First, there are no 
variables analogous to the p's and /^'s allowing a simple exact solution. The natural 
complementary variables are the Dirichlct data which are defined specially on each 
isospectral set rather than as a local function. The proper angle variables for this 
situation are discussed in terms of the theory of Abelian integrals on a suitable 
hyperelliptic surface 0[TOlll7l[23|24l|25l[3Ol|3l[46l|47]. 

Second, there is a sense in which the periodic situation is more subtle and in- 
teresting than the finite N structure. For finite N , the isospectral manifolds are 
open simplexes parametrized by {pj}^^ ('^'' P-j) with X]jLiPi=l- Topologically 
this is M^^^. For the periodic case, the isospectral manifolds are compact because 
n^=i '^3 rij^i Pj const, prevents aj or pj from going to zero (for OPRL one also 
needs boundedness of ^a|). By the general theory of completely integrable sys- 
tems, these isospectral sets are tori generically of dimension one-half the dimension 
of the symplectic manifold (p — 1 for OPRL, p for OPUC). 

The idea of the proofs of these theorems is to cut off the periodic matrix to a finite 
matrix and take limits of Poisson brackets. The limits of finite traces (normalized) 
will be the moments of the density of states, not the periodic eigenvalues — but we 
will see they are related. It will be useful to go back and forth between the two 
sets of fundamental symmetric functions: 

N 

E^' (11.19) 
E l<k<N (11.20) 

jl<---<ifc 

where we will need the following well-known combinatorial result: 
Proposition 11.3. For k < N, we have 

)(Ai, ...,Xn)^1 i'^^X,, ...,XN) + r (11.21) 

where r is a polynomial in {t^j^'^}^Zl and also a polynomial in {s'j^''}'^Zl . 

Proof. Here is a simple proof. Let A be the diagonal matrix with eigenvalues 
Ai,...,AAr. Let a'^(C^) be the fc-fold antisymmetric product and ®''(C^) the 
tensor product (see, e.g.. Appendix A of |39j). Let tt e E^, the symmetric group in 
fc-objects, act as on (g)'^(C^) via "Vyrixi lE" • ■ • ® Xk) — • • • (g) x^(fe). Then 



t)^'\X^,...,XN) = 

^k^\Xi, . . . , Xn) — 
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is the projection of ®*C^ to A*(C^). Thus, 

4'^)(Ai,...,Ajv) = Tr^HC«)(A'(^)) 



Trgj.(civ)(Pfc Cg)*^ A) 

^ ^ (-irTr(V,®'=yl) (11.22) 



fc! 

It is easy to see [SHI Appendix A] that if tt's disjoint cycle decomposition has 
one-cycles, ^2 two-cycles, etc., then 

Tr(V^ A) = Tr(A)^iTr(yl2)^2 . . . Tt[A''Y'' 

k 

= l[tf\x,,...,XNY^ (11.23) 

Since there are (fc — 1)! permutations with one fc-cycle, (lll.22p and (|11.23p imply 
pi.2ip with r a polynomial in t's. An easy induction then shows any tk is a 
polynomial in {sj}*^^]^ and so r is also a polynomial in the s's. □ 

The density of states of a periodic Jacobi or CMV matrix can be defined as the 
unique measure, d'^, (on M or dM) so that for all fee {0,1,2,... }, 



J Ja IP -^^ 



IT (11-24) 



Moreover, it is easy to see [IS] that if Jm{{o.k}^^i, {bk}k=i) is the cutoff Jacobi 
matrix obtained by taking the two-sided infinite Jacobi matrix and projecting onto 
the span of then for any fc G {0, 1, 2, . . . }, 

^ Tr{Jt) I X'^ d^{\) (11.25) 



2rn -f 1 

Similarly in the CMV case [H], if Cm is defined by taking q;_„i — > — 1 and — > /? 
(for any /3), then 

^ ■Tr(C,';J-^ / A'=d7(A) (11.26) 



2m + 1 

The key to the proof of Theorem 111.11 is the following. Define for {aj,f'j}j=i 
fixed, 

tu{e) = Tv{j{ef) (11.27) 

Sk{e)^ A,i(0)...A,,(0) fc = 0,l,...,p (11.28) 

jl<---<jk 

Theorem 11.4. Consider OPRL of period p. 

(i) For k < p ~ 1, Sk{0) is independent of 9. 

(ii) For k < p — 1, tk{9) is independent of 9. 

(iii) 

Sp{9) Sp{0) + {-l)p( f[ a A {2 -2 cos 9) (11.29) 
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(iv) 



tp{e) = tp{0) + {-l)Pp( Jl aj) (2 -2 cos 61) (11.30) 



In particular, for k <p, 

tk{0) = {-lY+\2p) Jl aj Skp + / A'= d7(A) (11.31) 

Remark. Directly from the form of the matrix (|11.3p . one can see (ii) easily, and 
(iv) with a little more work. 

Proof. By pT5)l and pT7)) . we see that 

p / ^ \ 

Y[ix-Xj{9)) = I J|aj j(A(a;) -2cos6') (11.32) 

which, expanding the products as 

f[{x - \,i9)) = + j2{-iy s,{9)xP'^ (11.33) 

immediately implies (i) and (iii). (ii) and (iv) then follow from (|11.2ip . and (|11.3ip 
then follows from (|11.24p . □ 

For OPUC we have: 

Theorem 11.5. Consider OPUC of even period p. 

(i) For fc = 1, 2, . . . , I — 1, 1 + 1, . . . ,p — 1, Sk(9) is independent of 9. 

(ii) For k — 1,2, ^ — 1, tk{9) is independent of 9. 
(iii) 

Sp/2{9)=Sp/2{0) + {-l)P/^(l[p,\2^2cos9) (11.34) 

^ 3=0 ' 

(iv) 

tviM = ip/2(0) + (-If/' f f n (2 - 2 cos 0) (11.35) 
/n particular, for k < p/2, 



3=0 



tki0) = {-ir/^+'p(f[p,)5,p/2+ fx'djiX) (11.36) 

^3=0 ^ '' 

Remark. It will suffice to have control of tj and ij for < \i\ < p/2 since Sp-j = Sj, 
as we will see. 

Proof. By pTH)) and pT7|) and the fact that zp/'^A{z) is a polynomial nonvanishing 
at z = 0, we have that 

f[{z - X,i9)) = P.) (^^/'A(^) - 2cos(0)zP/') (11.37) 

3 = 1 ^ 3=0 ^ 

which, via (|11.33p . implies (i) and (iii). (ii) and (iv) then follow from (|11.2ip . and 
pT36)) from ([1X24)) . □ 
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Proof of Theorem Let 

^(K,fe.}?=o) =Tr([J„({a.,6,}Lo)]''0 (H-SS) 

and 

Tt'^\WjMT=-^) = Tr([J({a,,6,},"L-™r) (11-39) 

So r^™^ is just T^™^ restricted to periodic sequences. 
Let 

tk{{aj,b,}P^,)^ I X'djiX) (n.40) 

As noted in (jll.25p . jm+i'^k"^'' ^ "^^^^ ^^^^ noting diagonal matrix 

elements of {Jm)^ are uniformly bounded in m for k fixed and equal to those of the 
infinite matrix so long as their index jo obeys |jo| < m — 2fc — 1. 
This same equality and polynomial nature show for 1 < j < p, 

.™r-^|- ("^«) 



^ y(m) 

doj+gp 



where {aj^"*, 5^"''} is periodic, and similarly for b derivatives. (|11.4ip holds so long 
as as |j + < m — 2fc — L It follows that 



,(0) 

6=6' 



a— a 

-hi") 



Since {f/™\f|™^} = 0, we conclude 



4 '^j 

{tfe,i,} = {n.43) 

for aU k,£. By pT3Tl) 

{tfe(0),i,(0)} = (n.44) 
for aU k, I. Thus by (|11.2ip . we see for 1 < A:, ^ < p, 

{sfc,s4 = (n.45) 

Since A(a;)-2 = (ai . . . a^)"^ 11^=1(^-^^(0)) and (|11.33|) holds and (ai . . . ap)^^ 
is in the Poisson center, we see that 

{A(a;),A(y)} = (n.46) 

Thus 

{A(x) - 2 cos(6'), A(y) - 2 cos(e'')} = (11.47) 

Evaluating this a,t x — Xj{9), y ~ Xk{d'), we get (jll.151) as usual. □ 

Proof of Theorem \11.2[ By repeating the arguments from the last proof, we see 
that (|11.44p holds and that 



{tfe(0),i,(0)} = {ife(0),i,(0)} = {UA 
From this and Theorem 1 11.5[ we see ioi 1 < k,£ < p/2, 



{sfe(O), s,(0)} = {sfe(O), s,(0)} = (11.49) 
By (jll.Sp . the product of the roots is 1, that is, 
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Since XjXj = 1 also, we see 

p 

so for 1 < fc, £ < p — 1, 

K(0),S£(0)} = (11.51) 
This, (|11.37p . and (fTTSll imply that 

|(n^^V(^)' (lip.VHj -0 (11.52) 

By a simple calculation if aj = jajje*^^, then 

so — log'dlj^o P'j) generates the Hamiltonian flow 

a J a^e'* (11.54) 

This transformation is implementable by a periodic unitary (see, e.g., |43] ) . and so 
it leaves the Xj{9) fixed. Thus, (|11.18[) holds, which implies 

iipv(\ipi)M^)]^o (11-55) 

j=0 ^j=o ^ ^ 
so (I11.52P imphes (|11.16p . which in turn implies pi.l7p . □ 

Remark. Unlike OPRL where Y[^=i '^j Poisson center, 11^=0 Pi generates 

a nontrivial flow, but one that leaves Xj{6) invariant. 

12. More Poisson Brackets for OPRL 

In this section, we will discuss some additional Poisson brackets for OPRL and 
related functions as a way of illuminating and extending the major results that 
we proved earlier. We want to begin by showing that one can go backwards 
from (|l.lip to (|1.15p / (|1.16p (or, more precisely, to the Poisson brackets (|1.16l) and 
{P„(x), P„(y)} = 0; that {Qn{x),Qn{y)} = is then a consequence of (|2.10p ). 

Clearly, 

{Pn{x), Pn{y)} = ^ {Xj,Xk} = (12.1) 

To see that 

{xj,pk} = i [6jkPj - PjPk] ^ (Ell) (12.2) 
we compute using p.5p and {xj,Xk} = 0, 

Pn{x) Pn{v) 



2{P„(x),Q„(y)} = 2^-{x„pfe} 



}J x~xj y~yk 

^ Pnjx) P,M ^ Pn{x)Pn{y) 

/ ,PjPk / .Pil T7 7 (12.3) 

^ ~ y~Xk ^ ' {x~ Xj){y -Xj) 
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The first sum gives Qn{x)Qn{y) by (I2.10|) and the second sum is {x 
y)"^[Pn{x)Q„{y) - Pn{y)Qn{x)] if we note that 



1 



1 



(x-xj) {y-xj) x~y 



(12.4) 



This estabhshes 

By (|2.1ip (with n replaced by n— 1) used in the first factor {Pn{x), P„(y)}, we get 
a relation among {P„_i(a;), P„(?/)}, {6i,P„(y)}, {al,Pn{y)}, and {Q„_i(a;), P„(?/)}. 
Since the first three are computed earlier, we get a formula for {(5„_i(x), P„(2/)}, 
namely, 



{Qn-l{x),Pn{y)} = 2P„_i(x) - 2Qn-l{x)~ 

_2/P„(a;)P„_i(2/) - Pn-i{x)Pn{y) 



x-y 



Pn-l{x)Pn-l{y) 



Next, we want to discuss Poisson brackets with 

Qn{x) 



mn{x) = - 



Pn{x) 



(12.5) 



(12.6) 



both to link to earlier work of Faybusovich-Gekhtman [13] and because one can 
then take n — > oo. 



Theorem 12.1. We have that 



{P„(x),m„(y)} = \ 
{Q„(a;),TO„(y)} = \ 



,{x)mn{y) 



Pn{x)mn{y) + Qn{x) 



-Qn{x)mn{y) +m„(j/) 



x~y 

Pn{x)mn{y) + Qn{x) 



x-y 



Proof. We have, since {P„(a;), P„(y)} — 0, that 

{P„(a;),m„(y)} = -P„(y)-i{P„(a;), Q„(y)} 
The second term on the right of (|1.16p leads to 

Pn{yy^Qn{x)Qniy) = ^Q7i{x)mn{y) 
The first term on the right leads to 

i Pn{yr'[Pn{x)Qn{y) " Qn{x)Pn{y)]{x - y)-' 

= (-^ Pn{x)mn{y) - i Qn{x)){x - y)"i 

proving (jl2.7p . 

On the other hand, since {Qn(a;), Qri(2/)} — 0, 

{Qn{x), niniy)} = Qn{y){Qn{x), -Pn(y)"^} 

^mn{y)Pn{y)-'{Pn{y),Qn{x)} 

= m„{y)P4y)-'{P„ix),Q„iy)} 
= mn{y){Pn{x),mn{y)} 



(12.7) 
(12.8) 

(12.9) 



(12.10) 
(12.11) 



so p^?71) implies ([T^ . (|12.10p follows from the symmetry of {Pn {x),Qn{y)} under 
x'^y, and (|12.1ip from □ 
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Theorem 12.2 ([13]). We have that 
{m„(z),m„(it;)} = i (m„(z)-m„(w)) 



TO„(z) - m„(w) 



z — w 



+m„{z)m„{w) 



(12.12) 



Proof. Since {P„(z), P„(u))} = = {Q„(z), (3„(w)} and {P„(z), Q„(w)} is sym- 
metric under z ^ w and {/,<?} is antisymmetric under f ^ g, 

Qn{z) 



{TO„(z),m„(ti;)} 



which is p2.12p . 



m„(z) 



TO„(z)m 

n 

(u;) 



{F„(z), (3„(w)} - {z ^ w) 
m,n{z) - mn{w)' 



(12.13) 

(z <-> Ul) 

□ 



As n ^ cxD, m„ has a hmit so long as {cj, 6j}j?:0 bounded (actually so long 

as the moment problem is determinate), so (|12.12p holds if m„(z) is replaced by 
m{z) for semi-infinite Jacobi matrices. 

The unnormalized transfer matrix has the form 

Pn Qn 

Of the sixteen Poisson brackets for these four functions, we have computed twelve. 
It would be interesting to know {Qn{z), Pn-i{w)}. Similarly, it would be interesting 
to know the Poisson brackets for the periodic transfer matrix and use them to prove 
(pTMl) . 

13. More Poisson Brackets for OPUC 

We begin by showing that one can go backwards from (|1.23p to (jl.lSp / p.l9p 
(or, more precisely, to (|1.19p and {Pn{z) , Pn{w)} — 0; {Qn{z),Qn{w)} then follows 
from the symmetry discussed after (|1.19p ). 

Clearly, 

{Pn{z),Pniw)} = ^ {ej,ek} = (13.1) 

To see that 

{0J,^J,k} = fJ-jSju - ^J-jfJ-k ^ <iL2M (13-2) 

we first need 



Lemma 13.1. For any distinct, z,w eC and e*^ G B, 



e'" + z e*" + w z + w / e'" + z 



e'" — z e'" — ui z — w \ e' 

Proof. Consider first e'^ = 1. Then 
1 + z 1 + w 



w 



(13.3) 



1 



2(z + 7«) 



and 



1 — z 1 — w 
1 + z 1 + w 



(l-z)(l-u;) 
2(z - w) 



1 — z I — w (1 — z)(l — w) 

which implies (|13.3p for e'^ = 1. 

In that formula, replace z by e~*^z and w by e~**w and thereby obtain the 
formula for general e'^. □ 
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Now we compute using ()4.9|) with zj — ^ , 

{Pn{z),Qn{w)} 

W + Zk 

z - z., - W - Zk 



Y^-{z„^,k}^Pni 



(13.4) 



i Z.J [5jklJ.j " MiMfe] ^"'^ 



w 



W + Zk 
W - Zk 









+ w\ 


1 gie. 






— w I 



= RHS of (fTTOl) 
verifying (|13.2p . In the above, we get (|13.5p by using 



- 1 



P„(z)P„H (13.5) 
(13.6) 

(13.7) 



and noting that the —1 term has j dependence only from [SjkiJ-j — fJ-jfik] which sums 
to zero for each fixed k. To get (|13.6p . we note that (|4.9p says 



Qn{z) 



J2 

k 



Pn{z) 



(13.8) 



and use Lemma 113.11 on the Sjk^JLj term. The first term in (|13.3p gives 
f^[P„(z)Q„(w) - g„(z)P„(w)] by (UMl). The second term in gives 
Pn{z)Pn{w) if we note that '^jkl^i^jk = ^ ^- IJ^j^-k term in (|13.5p 

gives the -Qn{z)Qn{w) by (|13.8p . 

Next, we consider Poisson brackets with Fn{z) given by (|4.8p and /n(z) — 

lr F„(^)-l n 

^ Lf„(z)+iJ- 

Theorem 13.2. We have 



{F„(z),F„(u;)} 



(P„(z)m„(w) + Qn{z)) 



-Pn{z) 

{Qn{z),F„{w)} = -Fniw){RHS of (UMl)) 
{F„(z),/„(u;)} = (l-u;/(7«))W^„(z,zi;) 

{Qn{z), fn{w)) = -{1 + wf{w))Wn{z, w) 



z + w 



i(z)m„(w) 



where 



Wn{z,w) = -- 



(Pn{z) + Qn{z)) - z{Pn(z) - Qn{z))fn{w) 



z — w 

Proof. (112. 9p is vahd with m„ replaced by F„. Thus 



(13.9) 

(13.10) 
(13.11) 
(13.12) 

(13.13) 



{P„(z),P„(u;)} = -- 



-Pn{wy 



z + w 



z — w 

+ Pn{z)Pn{w) - Qn{z)Qn{w) 

RHS of (fTXg)) 



(P„(^)Q«(U') -Pn(u^)gn(z)) 

(13.14) 
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(|12.11|) is also valid with rn„ replaced by Fn, so (|13.10p follows from (113. 9p . 
For the formula involving /„, we use (|4.28|) . Analogous to (112. lip , one finds 



and, by 



{C„(Z),/„H} = -wf{w){Sn{z)J{w)} 

{Sn{z), fniw)} = -W"^5„(?i))"^{S'„(z), C„(w)} 
' C„(z) + zSn{z)fn{w) ' 
Z — W 



(13.15) 

(13.16) 
(13.17) 
(13.18) 



= Wiw,z) 

by glSD/gH. 

We then get (fmTl) . (fmH using P„ = C„ + Sn, Qn = C„ - S'„ and (im5)l . 

(fms)) . □ 

Theorem 13.3 (Gekhtman-Nenciu 16J). FFe /laue that 



{F^{z),F,,{w)} = -- {F,,{z) - FAw)) 



[(r 


+ w\ 




— w ) 



{fn{z),fn{w)} 



(F^{z) - F„(w)) + 1 - F^{z)F^{w) 

f(z) - f{w) , f, . f, 

{zf[z) -wf{w)) 



(13.19) 
(13.20) 



z — w 

Remarks. 1. |16| has —i and —2i where we have — ^ and —i because their {•, •} is 
twice ours since, following [29j . they dropped the normalization of 37J which we 
keep. 

2. While we will separately derive (113. 19p and (|13.20p . it is an illuminating 
calculation to go from one to the other using F{z) = (1 + zf{z))/{\ — zf{z)). 

Proof. By the same calculation that led to (|12.13p . 

Q,i{z) 



{F^{z),F^{w)} 



(13.21) 









— w J 



(Fniz) - Fniw)) 

{z ^ w) 



+ l-Fniz)Fn{w) 

= RHs of imgi) 

where (|13.22l) follows from pTTO)) . 

Similarly, by (|4.28p and the same calculation that led to (|12.13p . 

iUz), Uw)} = -(zw)-' {5„(z), Cniw)} ~{z^w) 



(13.22) 



ifn{z) 



z — w 

i.fn{z) 



wSn{z)Sniw) 
Cnjz) £ Cn{w) 
_Sn{z) W Sn{w) 

2, 

Z.fn{z) Wfniw) 

W 



{z ^ w) 
— (z ^ w) 



(13.23) 



38 



M. J. CANTERO AND B. SIMON 

2 2^ J'^ ( ) 

{fn{z) - In{w)) -{z ^w) 

Z — W 

RHS of 113:201) 



□ 



Finally, we want to note that one can compute {$„(z), ^„(z)}, {$* (z), $*(?«)}, 
and {$„(z), {^'„(2), ^'^(m;)}. We did not know {$„(z), ^^(z)} or 

{\I'„(2;), $* (z)}. Nenciu [36], motivated by a preliminary version of this paper, 
has found this last missing bracket. With her formulae, one can get formulae for 
the brackets of the elements of the transfer matrix. 



Theorem 13.4. We have 

{<i>„(z), $„H} = {vi/„(z), ^,,{w)} = {K{z), KM} 

' z + w' 



{K{z),KM} 



{$„(z),*„H} 



Kiz)KH + Ki^)KH 

($„(z)^'„(u;) - *„(z)$„(u;)) 



(13.24) 



(13.25) 



(13.26) 



Proof. p.lSp and (|1.19p depend on (3 and are quadratic polynomials in (3. Equality 
for all /3 G 9D implies equality for all j3, and so equality of the coefhcients of the 
terms multiplying (3(3,(3, and 1. The (3(3 yield the three Poisson brackets 

and the terms with no (3 yield the three $„, Poisson brackets. □ 

Remarks. 1. The no (3 terms involve z$„(z) and w'^niw), but zw factors outs. 
The change of sign between (|13.25|1 and ()13.26p comes from the minus sign in ()4.1|) 
vs. the plus sign in ()4.2|) . 



2. It is an interesting exercise to use $^(z) = z" $„(l/z) to go from (|13.26p to 
(|13.25p . Since Poisson brackets of real functions are real, {f,g} = {/,<?}■ We get 
a minus sign from — z/2 = i/2, explaining the sign changes in the second and third 

terms. Because ( i/g-i/^ ) 



z+w ^j^g term does not change. 



Theorem 13.5. We have that 

{<i>n{z),KM} =iw 



$„(z)$;(w;)-$„(u;)<i>;(z) 



(13.27) 



z — w 

Remarks. 1. The same formula holds if $* is replaced by 5'„, ^* (since a„ ^ 
— a„ preserves {•,•}). 

2. The proof provides another proof of {$„(z), $„(u;)} = {$,* (z), $,* (w)} = 0. 

Proof. Define 

C„{z) = z$„(z) 5„(z) = $;(z) (13.28) 
Then Szego recursion becomes 

C„+i(z) = zC'„(2) - a„z5„(z) (13.29) 

5„+i(z) = 5„(z) - a„C'„(z) (13.30) 
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which has the same structure as (|4.16|) . (I4.17P except that ao, ai, . . . is replaced by 
ao, ai, . . . ,an- Moreover, at n = 0, 

Cn{z)wSn{w) — Cn{w)zSn{z) = ZW — ZW = 

So we can start the induction. Thus, by induction as in Section [5] {—i becomes 
the +i because of the complex conjugate change which flips signs of {an, an} to 
{an, an}), 

I \ Q I \\ - Cn{z)wSn{w) - Cn{w)zSn{z)' 

which is equivalent to (|13.27p . □ 

Using the recursion relations for <!>„ in terms of <I>„_i, one obtains 
Theorem 13.6. We have 

{$„_i(z), $„(w)} = ~ian-iw[ (13.31) 

\ z — w J 

. X ^, . U • ^ ^»-i(^)^;_i(w) - $.„,l(^;)$;_l(z) ^^ 

{$„ z w } = -izwl (13.32 

V z — w I 
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